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Abstract. We show that the singlet fraction ps and total magnetisation (or 
polarisation) m can bound the minimum concurrence in an ensemble of spins. 
We identify ps > (1 — m^)/2 as a sufficient and tight condition for bipartite 
entanglement. Our proof makes no assumptions about the state of the system or 
symmetry of the particles, and can therefore be used as a witness for spin entanglement 
between fermions. We discuss the implications for recent experiments in which spin 
correlations were observed, and the prospect to study entanglement dynamics in the 
demagnetisation of a cold Fermi gas. 
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1. Introduction 


Spin correlations have recently been observed in cold fermionic atoms as signatures 
of pairing, magnetism, and interaction strength [I112113]. It is interesting to ask if the 
observed correlations require pairwise entanglement. Since typical experimental samples 
contain thousands of atoms, full tomography is inaccessible; instead, one must find an 
entanglement witness based on a reduced set of measurements A commonly 

explored approach has been to measure spin squeezing p 0 El El [ini mi [12], however 
this approach is mainly limited to symmetric states or indistinguishable particles mw, 
and thus inapplicable to spin mixtures of fermions. 

An alternate characterisation may come from the degree of polarisation (or 
magnetisation) m, and the spin-singlet fraction ps of the ensemble (see Sec. [2] for precise 
definitions). The singlet state plays a key role in the physics of ultracold fermions, 
since a singlet spin wave function is required for s-wave interactions, which are the only 
interactions not suppressed at low energy by the centrifugal barrier. For spin mixtures 
near a Feshbach resonance [13], the pairing fraction can be measured by an adiabatic 
rapid passage that projects interacting pairs onto molecular dimers [Ml [151 [iH]- This 
enables direct measurement of ps in an ensemble. Singlet fraction is also proportional 
to the s-wave contact [T71 UHl [IHl [201 EH E21 EH EH ESI [2], and singlet pairs in an 
optical superlattice can also be mapped or projected onto excited motional states 

[261 EH EH EH [301 EH E]. 

It is well known that ps > 1/2 indicates pairwise entanglement in unpolarised 
(m = 0) ensembles [32] . The existence of a threshold is intuitive, since spin singlets are 
antisymmetric Bell states. Here we assume both m and Ps of an ensemble are measured, 
but make no assumptions about the form of the reduced two-body density operator 
which has 15 degrees of freedom. 

We hnd that the concurrence C of the ensemble can be bounded: 


C > max 


Ps - ^(1 -PsY - m2, 0 


( 1 ) 


This delineates a bound on the singlet fraction of an arbitrary two-body state that is a 
sufficient and tight condition for its entanglement, namely 
1 — 

This “singlet bound” is an extension of the Werner bound to partially polarised (m > 0) 
ensembles, where ps > 1/2 is sufficient but is not a tight bound. Our proof makes no 
assumptions about the state and therefore is a condition for bipartite entanglement in 
any ensemble of spins. Thus, Eq. ([H) can elucidate the nature of spin correlations in 
recent experiments with interacting ensembles of spin-half fermions, even in states far 
from equilibrium. 
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2. One- and two-body observables 

The spin state of pairs in a spin-1/2 ensemble can be described using the antisymmetric 
singlet state |so) = (Iti) “ lit))/V^ and the symmetric triplet states |to) = (|ti) + 
\i'\))/\/2, |ti) = |t_i) = These form an orthonormal set of basis states with 

well dehned angular momentum quantum numbers \S, S^)- The most general state can 
be written 

=Ps\so) {so\+ (giko)(ti|+gri^i)(so|)+ Y (3) 

iG{0,±l} ijG{0,±l} 

where the populations are normalised to Tr[p^^] = 1, and ps is the singlet fraction. 

The magnetisation is m = (rnx.,rny,mz) = Tr[Sp^^], where S = {a-^ + d'^)/2 
and are the usual Pauli spin operators. The reduced one-body states (e.g., 

Pa = TcbP^^) are completely dehned by a Bloch vector v. pa,b -^1/2 + va,b ■ d/2, in 
which I is the identity operator. Since m = \va + \vb-i 

= i i^A + ^1 + ‘^vaVb cos (3) (4) 

where jS is the angle between the two Bloch vectors, and m = |m|. 

In the problem we are considering, only the ensemble observables m and Ps are 
measured, not va, vb, or (3. One simple relation between m and ps is given by the 
normalisation of probability: 

Ps < 1 — m. (5) 

This can be shown by noting that = pu — P-i-i and Ps + YliPa — 1 ) from which 

the singlet fraction is bounded by Ps < 1 — |m^| — poo- Since \mz\ < m, Eq. ([5]) follows. 

3. Unentangled spins 

Let us start by Ending the singlet fraction of the separable state p^^ = p^ ® p^ where 
and p^ can be different mixed states. Since we are seeking a relation between two 
rotationally invariant quantities, m and Ps, we are free to choose the coordinate system, 
and align p^ along the ^ axis in Bloch space. Then p^ = p^t ITa) (TaI + PAi IJ-a) (J-aI, 
whereas p^ remains arbitrary, and we write it as p^ = Ns) Os I where i,j G 

{t, i}. The singlet fraction is 

Ps = \PA^C^ + \PAiCpf- ( 6 ) 

In terms of the Bloch vectors, p^t = (1 + s^)/2 and pa\. = (1 — '1^a)/ 2, whereas 
cx\. = {1 — Vb cos /5)/2 and c^tq = (1 + vb cos /9)/2, thus 

Ps = \{1 -vaVbcosI3). (7) 

Equation ([7]) has a simple interpretation for two pure states: when the Erst spin is along 
the +z axis of the Bloch sphere, the antiparallel (spin-down) fraction of the second spin 
is equally split between singlets and triplet zeros [33] . 
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For an ensemble of unentangled qubits, each of which is in the same unknown mixed 
state, Gisin noted that = (1 — m^)/4, and proposed measuring ps as a more efficient 
determination of m than measuring m [M]. We recover this result from Eq. ([7]) with 
va = vb and /3 = 0. However for arbitrary or Vb, we can only bound Ps'- eliminating 
/?, with Eq. (HD, 

p, = 1(1 - + \{v\ - 1 + - 1 )) < ^ ^ , ( 8 ) 


where the inequality holds because |nyi| < 1 and \vb\ < 1- Note that separable pure 
{va = ns = 1) states are examples of non-entangled states on the line = (1 — m^)/2, 
which demonstrates the tightness of Eq. (ED. (If however magnetisation is known only 
along z, but the full magnetisation possibly lies along another direction, the singlet 
bound Ps > (1 — nr^)/2 is sufficient but no longer tight.) 

We generalise the inequality in Eq. ([HD to all non-entangled states by considering 
a mixture of separable states i.e. PkPk^ where Pk is the probability of 

Pk^ — Pk ® Pk - The singlet fraction ps^ of each is still bounded by Eq. (|HD, thus 

1 - Efc Pk^l _ 1 - m2 


Ps ^ ^ PkPsk — 


(9) 


since m2 = J^kPk'i^l- Hence if the two-body state is non-entangled i.e. = 

Y.k PkPk ® Pk-: then Ps < (1 - m^)/2 holds. 


4. Concurrence of entangled states 

The contrapositive must also be true: if Ps > (1 — m^)/2, then is entangled. In 
fact, we hnd the concurrence [3H] of p"^'® can be bounded using ps and m, without any 
additional assumptions. 

First, we dehne a “spun state” as the state p^^ averaged uniformly over local 
rotations about the z axis, Uz{0) = U^{9) ® Uf{0)-. 

= deul(0)p-''>uM- ( 10 ) 

This transformation eliminates coherences between states in p^'® with different angular 
momentum quantum number Sz, since Uz{0) = exp[i6'S'^]. Populations and coherence 
between |so) and \to) are unaffected, leaving 

(p^^) =Ps\so){so\ + qo\so){to\+qo\to){so\+ ^ {pulU) {UD (11) 

iG{0,±l} 

which now has only six degrees of freedom. Crucially, because rotation can be 
implemented using local operation and classical communication (LOCC), the spun state 
is at most as entangled as the unspun state i.e. C{p^^) > C{{p^^)) [36] . 
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Figure 1. Each circle in singlet fraction ps vs. magnetisation m space is a randomly 
generated state described by Eq. (fTTl) . Blue squares have C = 0, while green 

circles have C > 0. Note that there are blue squares immediately beneath the singlet 
bound Eq. @ while there are none above, evidence that the bound is a tight and 
sufficient condition for entanglement. About the singlet bound, contour lines of C give 
the minimum concurrence (indicated by the colour scale) of a state with a given ps 
and m. For a given C, the line of minimum ps follows Eq. m from Ps = (1 +C)/2 at 
m = 0 to Ps = C at m = 1 —C. Several C values are given along physical limit (Eq. [5]). 


Next, we constrain the state to have polarisation m. Choosing the axis along the 
measnred direction of m, 


{p ) =Ps |so) (sol + a |to) (^ol + ko) (^o| 

I V / I ^ I V / I h — m , , , , 

+ ce \to) (sqI H --— |ti) (til H--— |t_i) (t_i| , 


( 12 ) 


where the normalised popnlations are Ps + a + b = 1 and the coherence is c = rjyjaps 
with 7] G [0,1]. 

Finally, we explicitly compnte the concnrrence of The eigenvalnes of the 

matrix R = {p^^) \/{p^^) i where = {cry®ay) {(Ty®ay) is the 

“spin-flipped” state, are 

1 r ,_ 1 1/2 

\i ^2 = —j= of + pI — 2c^ cos 20 ± y/ {a? + — 2c^ cos 20)^ — 4(c^ — apsY , (13) 

\/2 L -I 

A 3 = A 4 = |\/62 - nf. (14) 

The concnrrence is then 


C((p^^)) = max[0, Ai - A 2 - As - A 4 ] 

= max[0, \J{ps — ay + 4c^ sin^ 0 — \/6^ — m^], (15) 

which is nonzero when 

1 / 1 — 2a — A 

2 \ 1 — 2a -b 2ap^ sin^ 0/ 


( 16 ) 
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Since C{p^^) > C{{p^^)), Eq. ffTB]) provides a general bound on the singlet fraction for 
the entanglement of any with magnetisation m, triplet population a, and coherence 
c. With only the observables ps and m, this yields a sufficient condition for entanglement: 


Pi 


> sup 

a,ri,<j> 


1 - 2a 


m 


2 1 — 2 a + 2ap sin^ 


m 


(17) 


With a = 0 (which implies c = 0) and ps = (1 — m^)/2, we see that Eq. flT^ reduces 
to separable pure states, which fulfils Eq. ([ 8 ]) and saturates the bound. Another special 
case is the Werner state a = h/2 = (1 — Ps)/3, m = 0, and c = 0, for which Eq. (IT 6 |l 
becomes ps > 1 / 2 . 

The singlet bound found here (Eq. [2]) improves upon the generalised witness of 
Ref. [2Z], which when applied to (p^^) with a = 0 , yields the sufficient condition 
Ps > (1 +V 1 + 3C2)/3. 

The bound can be generalised to a threshold for finite concurrence, knowing only 
m and ps, by noting that the minimum of Eq. flT^ occurs when a = 0. Along with the 
constraint of a normalised probability, ps + b = 1, this gives Eq. o. Solving for ps, this 
gives a tight and sufficient condition for p^^ having at least concurrence C, namely 

1 _ (J 2 _ ^2 

Ps >- 7 — and m < 1 — C (18) 

^ 2(1-C) - ^ ^ 

where Eq. ([2]) is now found from the condition C > 0. Equations ([T]) and m are the 
central results of our work. 

We verify these relations by generating random mixed states that span the Ps and 
m space, and computing their concurrence. Each point in Fig. [1] corresponds to one of 
five thousand random spun mixed states. The blue squares have C({p^^)) = 0 and are 
not entangled while green circles have C({p^^)) > 0 and are entangled. All points lie 
within the physical limit, Eq. ([5]). The absence of non-entangled states above the singlet 
bound demonstrates that Eq. ([2]) is a sufficient condition for entanglement of while 
the existence of non-entangled states immediately beneath the bound demonstrates the 
tightness of the condition. Note that there are also entangled states below the singlet 
bound, as it is not a necessary condition for entanglement. Figure [T] also shows contour 
lines of minimum C determined from a larger set of random matrices. The locus of 
points with at least concurrence C or greater is bounded by Eqs. flTSl) and ([5]). 


5. Discussion 

Several recent experimental works can be re-interpreted in light of our results. We 
will consider three measurements sensitive to p^; mapping onto vibrational states in a 
superlattice, sweep-projection onto singlet dimers, and measuring the s-wave contact. 
We focus on experiments with fermions, even though our results apply to mixtures with 
any exchange statistics. 

Controlled collisions in optical superlattices have been used both to create and 
to detect pairwise entanglement [261 EH EHl EHl EHl EH H]. However, when the effect of 
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uncontrolled collisions are measured with the same technique, the efficiency of observing 
Ps may be hampered by a randomised choice of pairs, if a simple lattice is pairwise 
projected into the superlattice. For instance, Greif et al. [1] hnd that in a dimerised 
lattice, the singlet fraction of fermion pairs is at least Ps = 0.31. This was an effective 
probe of spin correlations, but insufficient to prove entanglement by Eq. fl2l) . 

The association of atomic fermions into s-wave pairs can also requires an initially 
singlet spin state. Thus efficiency V of association is a lower bound on pg. For example, 
sweeping the magnetic field across a Feshbach resonance in experiments with unpolarised 
Fermi gases of and ®Li, V as high as 85% is observed [38l [Ml US] . This surpasses the 
50% upper limit discussed in Refs. [IHlIlO] which is also seen as an apparent limit in some 
experiments [Hill]. We interpret this limit as Ps = 0.5, which is the maximum singlet 
fraction of a non-entangled state: experiments (and theoretical treatments) Ending 
V < 0.5 use separable states, whereas experiments observing V > 0.5 allow multiple 
collisions to occur before or during the magnetic field ramp. In some conditions, these 
collisions have produced pairwise entanglement. From Eq. ([T]), we can infer that the 
concurrence was C > 0.7 for Ps ^ 0.85 in Refs. [39l EHl US]. 

Pairwise-entangled states of an unpolarised Fermi gas are not surprising: in a weakly 
interacting Fermi s-wave superfiuid, each spin-up fermion is (monogamously) entangled 
with a spin-down partner. However entanglement dynamics in a polarised gas is an 
active topic of discussion. Calculations of the s-wave contact X in a polarised Fermi 
gas [21 02] have shown that X oc 1 — at high temperature, and X cx 1 — m at low 
temperature. Since X reflects interaction strength, which in turn requires spin-singlet 
wave functions between fermions, this is similar to a study of Ps versus m. The conversion 
of X to an absolute value of ps requires many-body theory and precise knowledge of 
density, temperature, and interaction strength. For this reason the spin correlations 
found by Bardon et al. [2| using X and m, for instance, cannot easily be classified 
using the singlet bound. More clear would be to study demagnetisation dynamics using 
association efficiency V. One would anticipate a temperature threshold, below which 
the gas evolves from a separable state to a pairwise-entangled state through random 
collisions. 

In sum, we have established a sufficient and tight condition for bipartite 
entanglement between spin degrees of freedom in an arbitrary system of spins, without 
any assumption of equilibrium, population balance, or symmetry. We find that the 
concurrence can be bounded simply by the magnetisation and singlet fraction, through 
Eq. ([1]). This enables the distinction between classical spin correlations and necessarily 
quantum correlations in ensembles of ultracold fermions. 
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